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Abstract 

A summary of measurements of the fragmentation of charm quarks into a specihc 
hadron is given. Measurements performed in photoproduction and deep inelas¬ 
tic scattering in e^p, pp and collisions are compared, using up-to-date 

branching ratios. Within uncertainties, all measurements agree, supporting the 
hypothesis that fragmentation is independent of the specihc production process. 
Averages of the fragmentation fractions over all measurements are presented. 
The average has signihcantly reduced uncertainties compared to individual mea¬ 
surements. 
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1 Introduction 


The production of specific charm hadrons has been measured in different regimes and en¬ 
vironments: in e’''e“ collisions at i?-factories CHS] and in Z decays [9HI3], in collisions 
in photoproduction (PHP) [IH[IS], deep inelastic scattering (DIS) [IBHIB] and in pp co¬ 
llisions [MM]. 

The fragmentation process is soft and hence can not be calculated with the techniques 
of perturbative QCD (pQCD). Therefore, these measurements are a necessary ingredient 
for any QCD prediction of charm-hadron production. In this context, it is important to 
validate the hypothesis that fragmentation fractions are universal, i.e. independent of the 
hard production mechanism. Thus, once precisely measured in one experiment, they can 
be applied in any reaction. Another important check is that the snm of fragmentation 
fractions of all known weakly decaying charm hadrons is equal to unity, thus checking if 
all weakly decaying states are known. 

To achieve these goals, a comparison of fragmentation-fraction measurements obtained 
in different prodnction regimes is performed using a combination of individnal measure¬ 
ments. Due to independent data sets and different detector types and constructions, 
the experimental statistical and systematic uncertainties in most cases can be treated 
as uncorrelated between measurements. However, a careful treatment of correlated nn- 
certainties due to common usage of branching ratios and theory inputs is essential, as 
in many measurements these are one of the leading uncertainty sources. In the past 
several combinations of fragmentation-fraction data were performed with fewer inputs: 
the summary of the charm fragmentation fractions in e'''e“ at the Z resonance iza, the 
combination of e’''e“ measurements [261 I^T] as well as the combination of and e^p 
measurements [2S]. Compared to those, the present analysis extends to a larger set of 
measurements, in particnlar the hnal measurement in PHP by the ZEUS experiment at 
HERA [151 , the pp measurements from LHCb [MI2U] . ALICE [2IH23|, ATLAS [21| and the 
A+ measurements from the BABAR experiment [8|. It uses the up-to-date branching-ratio 
valnes [MISHUSO], treats correlations of branching-ratio uncertainties and recent theory 
predictions with reduced uncertainties [211132] as input. 


2 Combination of individual measurements 

2.1 Update of input measurements to recent branching ratios 

To make separate inputs consistent, the original measnrements are corrected to the same 
up-to-date world averages of branching ratios of the charm-hadron decays, H, summarised 
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in Tab. [H Most of the values were taken from Ref. [25]. The B{D*^ —)■ was 

calculated from the two most precise measurements [2^130] of B{D*^ — /B{D*^ —)■ 
L)° 7 ) assuming B{D*^ + B{D*°^ D^-f) = 1. 


2.2 Calculation of the fragmentation fractions 

In this paper the charm-quark fragmentation fraction to a specific hadron H is defined as 
the production cross-section via charm quark divided by the production cross-section of 
the charm quark: 

f{c^H) = a{H)/a{c). (1) 

The charm hadrons produced in the decays of beauty hadrons are not considered. The 
Standard Model makes precise predictions for the total charm cross-section in e’''e“ co¬ 
llisions, therefore, for those processes it is possible to calculate /(c —)■ if) according to 
Eq. dl]). Sufficiently precise predictions for the charm-quark production in pp and e^p 
collisions are not available. However, it is possible to make an assumption that the sum of 
charm-quark fragmentation fractions to all known weakly decaying charm hadrons in the 
end of the fragmentation process is unity. Then the charm-quark fragmentation fraction 
to a specific hadron can be calculated as the ratio of the hadron-production cross-section 
over the sum of cross-sections of all known weakly decaying (w.d.) charm hadrons 

nc^H) = a(H)/J2a(H). ( 2 ) 

W.d. 


To obtain the charm-quark fragmentation fractions according to Eq. (]2|), in addition to the 
production cross-sections of D mesons and A+, it is necessary to know the cross-sections 
of the weakly decaying and states. Those states are poorly studied, therefore as in 
Ref. [H] it is assumed that ratios of fragmentation fractions of charm and strange quarks 
into the corresponding baryons are similar, /(c H+)//(c A+) = /(c —)■ H°)//(c —)■ 

AJ) = f(s ^ H-)//(s ^ A») and /(c ^ !d2)//(c ^ AJ) = f(s ^ !J-)//(s ^ A»). In 
this approach the sum of the production cross-sections of these states can be estimated 
as 

+ cr(H°)+ a(f2°) = Act(A+), (3) 

where we define 


a ^ /(s -»• 5 ) /(g -»■ n ) 

/(g ^ AO) ^ /(g ^ AO) 


0.136 ±0.006. 


(4) 


The value of A is calculated using the most precise set of s quark fragmentation fractions 
/(g ^ “-) = 0.0016 ± 0.0003, /(g ^ 12-) = 0.0258 ± 0.0010 and /(g ^ AO) = 0.3915 ± 
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Decay 

In this work (%) 

In experiments (%) 

D*+ D+ 

32.30 ±0.50 [25] 


D*+ D+7r° 

30.70 ±0.50 [25] 

30.70 ±0.50,BELLE [7] 

D*+ D°tt+ 

67.70 ±0.50 [25] 

68.13 ± 1.40,ALEPH [IT] 67.70 ± 0.50,ALICE I2TH231 
55.00 ± 4.00,ARGUS [SHS] 67.70 ± 0.50,ATLAS [21] 
67.70 ± 0.50,BELLE [7] 67.60 ± 0.50,H1 [IS] 

67.70 ± 0.50.ZEUS I14H161I33I 67.70 ± 0.50,LHCb l20l 

68.30 ± 1.40,OPAL PITOl 

D*0 jjO 

100.00 I29II30I 

100.00,H1 [18] 100.00,ZEUS [16] 

D*0 D°tt° 

64.94 ± 0.89 12911301 

61.90 ± 2.90,BELLE [7] 

Dt+ ^ Z)+7 

93.50 ±0.70 [25] 

94.20 ± 0.70,BABAR [6]|8] 

D+ K-1T+TT+ 

9.46 ± 0.24 [25] 

9.13 ±0.19,ALICE [2TH23] 7.70 ± 1.00,ARGUS [SHS] 
9.13 ±0.19,ATLAS [24] 9.20 ± 0.60,BELLE [7] 

9.00±0.60.H1 [THl 9.51±0.34.ZEUS I16II33I 

9.20±0.60,ZEUS [M] 9.13 ± 0.19,ZEUS [B] 

9.13 ± 0.19,LHCb 179112(11 

D+ K°Tr+ 

1.53 ±0.06 [25] 


D° K-TT+ 

3.93 ± 0.04 [25] 

3.85 ± 0.09,ALEPH [77] 3.87 ± 0.05,ALICE 127112^ 
3.71 ± 0.25,ARGUS [SHS] 3.88 ± 0.05,ATLAS [24] 
3.80± 0.09,BELLE [7] 3.83 ± 0.09,H1 [18] 

3.80 ± 0.07.ZEUS I16II33I 3.80 ± 0.90.ZEUS 1711 

3.88 ± 0.05,ZEUS 17^ 3.89 ± 0.05,LHCb I191I20I 

—>■ i^“7r"*"7r"*"7r~ 

8.07 ±0.23 [25] 


D+ ^ K*°{K-tt+)K+ 

2.61 ±0.09 [25] 


Dt <j){K+K-)-K+ 

2.27 ±0.08 [25] 

2.28±0.12,ALICE 121 11221 2.16 ± 0.28.ZEUS ITBl 

2.28 ±0.12,ZEUS 1751 0.00,LHCb I191I20I 

D+ </)7r+ 

4.50 ± 0.40 [25] 

3.60 ± 0.90,ALEPH [77] 3.60 ± 0.90,BELLE [7] 
3.60±0.90,H1 [78] 3.60 ±0.90,ZEUS [H] 

Dt K+K--K+ ,7Me\ 

1.83 ±0.06 I251I34I 

1.85 ±0.11,ATLAS [24] 

Dt K+R'-7r+,20MeV 

2.23 ± 0.08 I251I34I 

2.24 ±0.13,LHCb I191I20I 

—>■ 7r"^7r“ 

69.20 ±0.05 [25] 


K>0 K+TT- 

66.67 [25] 

66.67,DELPHI [77] 

(f) K+K- 

48.90 ±0.50 [25] 

49.10 ± 0.80,ALEPH [77] 49.00 ± 1.00,ARGUS [3}[5] 
48.90 ± 0.50,BABAR [6]|8] 49.10 ± 0.60,BELLE [7] 
49.00 ± 1.00,CLEO P[2] 49.10 ± 0.80,DELPHI [77] 
49.20 ± 0.70,H1 [78] 49.10 ± 0.60,ZEUS 1731 

49.10 ±0.80,OPAL 111117171 

—>■ P'K'^ 

63.90 ±0.50 [25] 


A+ pK--K+ 

6.841°;“ [25] 

5.00 ± 1.30,BABAR [6]|8] 5.00 ± 1.30,BELLE [7] 

5.00 ± 1.30,ZEUS I141I1SII17I 5.00 ± 1.30,LHCb [75] 

AJ —>■ A^(p7r"*")7r''" 

0.95 ±0.10 ±0.05 [25] 

0.68 ±0.18,ZEUS [77] 

A^ —)■ A®7r"^7r"*"7r“ 

3.59 ±0.28 [25] 


^ A°7r+ 

1.46 ±0.13 [25] 


A+ E.-K+TT+ 

0.70 ±0.08 [25] 


A+ ^ pK° 

3.21 ±0.30 [25] 


A^ pK^{fK^'K~) 

1.09 ±0.06 ±0.11 [25] 

0.80 ±0.21,ZEUS [IT] 

A^ —)■ pK^'K^'K~ 

3.50 ± 0.40 [25] 



Table 1: Branching ratios used for calculations. The second uncertainty for the 
B{Af) is the uncertainty of decay branching ratios of daughters. The numbers in 
the Df decay branching ratio indicate the used \M{K^K~^) — M((/)(1020))| mass 
windows. For the experiments which measured combination of cross-sections and 
branching ratios, the values of branching ratios are not given. 














































0.0065 from Ref. [25] obtained at LEP. Hereby, the sum of production cross-sections of 
all weakly decaying states is 

= a{D°) + ct(R>+) + ct(R>+) + a(A+) + Aa(A+). (5) 

w.d. 

The fragmentation fractions calculated according to Eq. ([T]) for the collisions and Z 
decays allow an independent check that 

S = /(c ^ £>») + f(c ^ B+) + /(c ^ B.+) + /(c ^ A^) + A/(c ^ A+) (6) 
is close to unity with sufficient accuracy. 

2.3 Combination procedure 

The combination of the measurements used in the present analysis is based on numerical 
minimisation with respect to observables of interest. The numerical minimisation was 
performed with the MINUIT package [35] and the procedure for calculation of itself is 
outlined below. 

For a set of m measurements and corresponding expectation values calculated from n 
parameters, a column-vector of the residuals i?(l x m) is calculated as a difference of a 
measurement and the corresponding expectation. The covariance matrix, l/(m x m), is 
calculated as 

Vij = 5ij + ^ ^ CjxCi^k, 

k 

where Ui stands for an uncorrelated uncertainty of i-th residual, Ci^k stands for the corre¬ 
lated uncertainty of source k of the Tth measurement and the sum runs over all sources 
of correlated uncertainties. The is then calculated as 

= R^V-^R. 

The correlated uncertainties are treated multiplicatively in the construction of the covari¬ 
ance matrix, i.e. the relative uncertainties are used to scale the corresponding expectation 
values instead of the measurement. This avoids the bias for normalisation uncertain¬ 
ties, such as branching ratio uncertainties, which are the main correlated uncertainties 
considered in the presented analysis. The statistical and uncorrelated systematic uncer¬ 
tainties are treated additively. Data sets and their systematic uncertainties are assumed 
to be independent between experiments. In addition, most of the measurements do not 
contain the information about a potential correlation between cross-section values for 
different charm hadrons. Therefore, in the following all experimental uncertainties are 
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treated as uncorrelated, unless otherwise stated. Uncertainties on the combined values of 
the fragmentation fractions are determined using the Hessian method with the criterion 

= 10 . 

The evaluated total uncertainties on the free parameters comprise experimental, branching 
ratio and uncertainties of the A parameter. 

The combination of all the measurements is obtained imposing the normalisation con¬ 
straint on the sum of all ground state hadrons by adding an additional “measurement” 
of S calculated from Eq. ([6]) with an uncertainty on A and the corresponding prediction 
S = 1. In order to keep the main result with the normalisation constraint as model 
independent as possible, no theory inputs on the charm cross-section are used in such a 
combination, and the fragmentation fractions are calculated according to Eq. ([5]). For 
the same reason, any measurements that require theoretical inputs for conversion into 
cross-sections or fragmentation fractions and do not have these inputs in the original pub¬ 
lications are also excluded from the main combination. However, such data are included 
in a more constrained combination. In the following, treatment of such measurements 
will be discussed case-by-case in the relevant sections. 

The quantities commonly used as Monte Carlo generator parameters, 

_ /(c —)• cu) ^ /(c —»• jA°) — f{c ^ D*^)Bd*+^do 
“ f{c ^ cd) ~ f{c ^ D+) + f{c ^ ’ 

^ 2/(c^ cs) ^ 2 /(c^T>+) 

/(c^ch/dV f{c^D+) + f{c^D^y 

^ 2f{c^cs) = ll ~ 2/(c^ T>;+) 

f{c ^ cu/d) ^ ^ ^ /(c ^ D*+) + f{c ^ D*0) 

and 

, ^ /(c ^ cu/d){J = 1) /(c ^ D-*) + /(c ^ £>•") 

'' f(c^cu/d)( J = 0)~ /(c ^ £>+) +/(c ^ 

were calculated from the £t results with the full error propagation and taking into account 
the correlation between parameters. 

^ As an illustrative example, when for a given combination set-up the inputs are the m cross-section 
measurements. These data define the mxm covariance matrix V with the experimental statistical and 
systematic uncertainties contributing to the diagonal elements and the correlated uncertainties setting 
the off-diagonal elements and contributing to the total uncertainties on the diagonal. The correlated 
uncertainties considered are those related to A in Eq. m and branching ratios. The residuals are 
obtained subtracting from the measurements cross-section expectations calculated from n free param¬ 
eters in the fit, which could be fragmentation fractions, total charm cross-sections, kinematic factors, 
etc. The details of this calculation are outlined in the each section. With all these components at hand 
the can be evaluated and iteratively numerically minimised with respect to the free parameters. 
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3 Charm-quark fragmentation into hadrons in e^e 
collisions 

Measurements of charm-hadron-production cross-sections in e^e“ collisions hadrons were 
based on the differential momentum spectrum dcr(e’''e“ H)/dxp. The extrapolation to 
the total cross-section was made in the orimnal papers using a theoretical fragmentation 
function (e.g. Bowler [SB] or Peterson [27]) Q 

As mentioned before, the precise predictions of the total charm-production cross-section 
in allow calculation of the fragmentation fractions without constraints on the sum 

of fractions. This way the used hypothesis about the sum of fragmentation fractions 
(Eq. (]6|)) can be verihed. 

3.1 Charm-quark fragmentation fractions from measurements 
at 5-factories 

The i?-factories provided many results on charm-hadron production around the T reso¬ 
nances, which can be used for the calculation of the charm-quark fragmentation fractions 
in hadrons (see Tab. [2|). The results of the CLEO [Tl|2] and ARGUS [3H5] experiments are 
represented as a product of the charm-hadron cross-sections times decay branching ratios, 
a{e'^e~ —?■ H) ■B{H —)■ daughters). The BELLE experiment [7] provided measurements of 
H). The BABAR experiment [S] provided a measurement of an average number 
of A+ —)■ pK~7i^ decays per hadronic event 

NZ . B(A- ^ PK-.*) = X B(At pK-.*) = 

= 2R, ■ f{c ^ A+) • S(A+ ^ pK-TT^), 

where Rc = average number of charm-quark pairs per hadronic event. 

A prediction of Rc is needed to use this measurement as an input, therefore, as discussed 
in Sec. 12.31 it is used only in the case of a{e^e~ —)■ cc) hxed to a theoretical prediction. 

^ A proper extrapolation procedure requires only the hadrons produced directly in fragmentation to 
be used in the fits. The hadrons produced in decays of excited charm hadrons should be treated 
separately. 

In many cases the limited precision of the measurements makes this requirement hard to follow and the 
decay part of the meson production is treated together with the fragmentation part. In the cases where 
the contribution of hadrons from decays is comparable with the contribution of the direct production 
in fragmentation, e.g. for and the joint treatment could bias the results. 
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Decay 

CLEO 

(T ■ B, pb 

ARGUS 

cr, pb 

BABAR 

cr ■ B, pb 

BELLE 

cr, pb 

D*+ D°it+, 





D° K--K+ 

17.0 ±1.5 ±1.4 [3 

690 ± 80 [3] 


598 ± 2 ± 77 [3 

D*+ D+tt°, 





D+ K--K+-K+ 




590 ± 5 ± 78 [3 

D*+ D°-k, 





—>■ /C“7r“7r"’"7r"*" 

33.0 ± 3.0 ± 1.8 [T] 





30.0 ± 3.0 ± 5.8 [3 




D*o ^ £)0^o 




510±3±84[3 

D° K-1T+ 

52 ± 5 ± 4 [3 

1180 ± 150 [3] 


1449 ± 2 ± 64 [3 

D+ K~n+n+ 

51±7±2[3 

650 ± 90 [3] 


654 ± 1 ± 36 [3 


(7 ■ B, pb 

(T ■ B, pb 

N'g ■ B X 10'^ 

cr, pb 

A+ ^ pK-TT+ 

13.5 ± 4.0 ± 1.4 [3 

9.0 ±1.2 ±1.0 [5] 

284 ± 4 ± 9 [8] 

189 ± 1 ± 66 [3 

—>■ pK~7v^ 

10.0 ±1.5 ±1.5 [2] 




A+ ^ K°p 

4.6 ± 0.6 ± 0.8 [2] 




—>■ K^p'n-^7r~ 

4.5 ±1.3 ±1.1 [2] 




A^ A°7r+ 

1.9 ±0.3 ±0.3 [2] 




A^ —>■ A*^7r"*"7r"*'7r“ 

6.8 ± 1.0 ± 1.3 [2] 




A+ ^ E.-K+TT+ 

1.6 ±0.4 ±0.3 [2] 





(T ■ B, pb 

(T ■ B, pb 

cr ■ B, pb 

cr, pb 

D+ 4>Tr+ 

7.2 ± 1.9 ± 1.0 [3 

7.5 ± 0.8 ± 0.7 g] 

7.55 ±0.20 ±0.34 [6] 

231 ± 2 ± 92 [3 

D* —> (/)7r+7 



5.8 ±0.7 ±0.5 [6] 



Table 2: Measured cross-sections, a, cross-sections times branching ratios, a ■ 
B, and N'^ ■ B(see text for explanation) for the production of the charm hadrons 
produced in e'^e~ collisions at centre-of-mass energies of y/s = 10.55 GeV fJ]/, s/s = 
10.52 GeV 00/, = 10.47 GeV ^/s = 10.5 GeV ij, y/s = 10.2 GeV and 

y/s = 10.54 GeV The numbers are given as in the original publications with 

uncertainties related to branching ratios omitted. The first uncertainty is statistical 
and the second is systematic. 
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For the calculation of the charm-quark fragmentation fractions a ht procedure is used 
as described in Sec. [2j The total charm-quark-production cross-section is calculated as 
described in App. The fit parameters are the fragmentation fractions. The obtained re¬ 
sults are given in the middle column of Tab. [3l The sum of the charm-quark fragmentation 
fractions into weakly decaying states calculated according to Eq. (|S]), S't = 0.9701 ±0.0284, 
is consistent with unity. 



Fixed a{e^e —)• c) 

Constrained S 

f{c^D*+) 

0.2470 ± 0.0137 

0.2525 ±0.0155 

f{c^D*0) 

0.2241 ± 0.0304 

0.2291 ±0.0316 

f{c ^ D*+) 

0.0532 ± 0.0082 

0.0544 ± 0.0085 

/(c^D+) 

0.2639 ±0.0139 

0.2698 ±0.0125 

f{c ^ D°) 

0.5772 ± 0.0241 

0.5901 ±0.0140 

f{c^D+) 

0.0691 ± 0.0045 

0.0707 ± 0.0048 


0.0526 ± 0.0031 

0.0611 ±0.0060 


19.2 

17.0 

^dof 

21 

20 

s 

0.9701 ± 0.0284 

1.0000 ±0.0005 

^u/d 

0.9508 ± 0.0752 

0.9508 ± 0.0752 

pd 

0.5601 ± 0.0432 

0.5601 ± 0.0431 

7s 

0.1644 ±0.0121 

0.1644 ±0.0121 

7l 

0.2257 ± 0.0385 

0.2257 ±0.0385 


Table 3: Average of charm-quark fragmentation fractions in hadrons in e'^e~ 
collisions around y/s = 10.5 GeV. The quantities S, Ru/d, Py Is were recal¬ 
culated from the fit results taking into account correlation of fit parameters. The 
value of minimised o,nd the number degrees of freedom of the fit Udof are given 
as well. 


The combination is also done according to Eq. ([2]) and imposing the constraint At —1 = 0, 
to be consistent with the definition used for e^p and pp data. The fit parameters 
are the fragmentation fractions and the total charm cross-section. The centre-of-mass 
energy dependence of the charm-quark cross-section is accounted for, according to for¬ 
mulae in App. 1^ taking the total charm-quark cross-section at a centre-of-mass energy 
^/s = 10.5 GeV as a reference. The results are given in Tab. [3] (right column). In this 
approach, the precise BABAR measurement of A^ production [S] is not included in the 
combination since it requires usage of the Rc theoretical prediction. The latter has an 
influence on other fragmentation-fraction results. 






















3.2 Charm-quark fragmentation fractions from measurements 
at LEP 

The LEP collider provided many results on the charm-hadron production. The most 
valuable for the studies of fragmentation are results obtained from Z decays. Most of those 
results are represented in the form of fraction of charm events multiplied by branching 
ratios r(J.^hadrons) -^) ■ daughters) (see Tab. 0]). In addition, ALEPH pT] . 

DELPHI 113] and OPAL (TU] provided measurements of /(c —?■ D*+) from the hts of 


Particle or decay 

OPAL 

(10-6) 

ALEPH 

(10-6) 

DELPHI 

(10-6) 

D+ <I>{K+K-)it+ 
D+ K*°{Kn)K+ 

560 ± 150 ± 70 [9] 

352 + 57 + 21 nn 

765 ± 69 ± 37 [12] 

624 ± 122 ± 73 [12] 

D° K-TT+ 

3890 ± 27012 “ 0 

3700 ± 110 ± 230 [lU 

3570 ± 100 ± 146 [H] 

D+ K-iT+n+ 

3580 ± 46 OI 310 [9] 

3680 ± 120 ± 200 [U] 

3494 ± 116 ± 140 [H] 

A+ pK-TT+ 

410 ± 190 ± 70 [9] 

673 + 70 + 37 nn 

743 ± 155 ± 78 [I2] 

D*+ D°{K-n+)n+ 

1041 + 20 + 40 ITOl 


1089 + 27 + 39 [T 2 I 


/(c ^ H) 

(10-3) 

f{c ^ H) 

(10-3) 

f{c ^ H) 

(10-3) 

D*+ 

222 ± 14 ± 14 [TO] 

233.3 ± 10.2 ± 8.41 [TT] 



f{c 

(10-3) 

f{c 

(10-3) 

f{c 

(10-3) 

D*+ D°-k+ 



174 ± 10 ± 4.2 [13] 

D*+ <p[K+K-)-K+'y 


69 ± 18 ± 7 [n] 



Table 4: LEP measurements of the products of the partial decay width of the 
Z into cc quark pairs, charm-hadron-production fractions, f{c —)■ H), and 

corresponding branching ratios. The numbers are given as in the original publica¬ 
tions with uncertainties related to branching ratios omitted. The first uncertainty 
is statistical and the second is systematical. 

fragmentation functions (see Tab. 0]). 

For the calculation of charm-quark fragmentation fractions, a ht procedure is used, as de¬ 
scribed in Sec. [21 The theoretically calculated value that is used, = 0.17223 ± 

0.00001 [32], is in agreement with the experimental world average 0.1721 ± 0.003 [25]. 
The ht parameters are the fragmentation fractions. The results are given in the middle 
column of Tab. |5l The sum of the charm-quark fragmentation fractions into weakly de¬ 
caying states calculated according to Eq. ([H|), 5’^ = 0.9292 ± 0.0261, differs from unity by 
2.7 standard deviations. 

The combination is also done using Eq. ([2]), and imposing the constraint Sz — I = 0, 
to be consistent with the dehnition used for e^p and pp data. The ht parameters are 
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Fixed 

i had 

Constrained S 

f[c^D*+) 

0.2369 ± 0.0064 

0.2454 ± 0.0071 

f{c^ D*+) 

0.0545 ± 0.0144 

0.0547 ± 0.0145 

f{c^D+) 

0.2267 ±0.0100 

0.2429 ± 0.0102 

f{c^ DO) 

0.5470 ± 0.0215 

0.5894 ± 0.0132 

f(c^D+) 

0.0925 ± 0.0082 

0.0996 ± 0.0083 

/(c^ A+) 

0.0555 ± 0.0065 

0.0600 ± 0.0066 


6.7 

7.8 

^dof 

13 

13 

s 

0.9292 ± 0.0261 

1.0000 ± 0.0005 

^u/d 

0.9987 ± 0.0627 

1.0348 ± 0.0580 

pd 

0.6119 ±0.0185 

0.6000 ± 0.0177 

7s 

0.2390 ± 0.0224 

0.2394 ± 0.0223 


Table 5: Average of charm-quark fragmentation fractions into hadrons in Z 
decays. The quantities S, Ru/d, Py 7s recalculated from the fit results 
taking into account correlation of fit parameters. The value of minimised and 
the number degrees of freedom of the fit Udof are given as well. 

the fragmentation fractions and the r(z-^hadrons) ^^^io. The results, given in Tab. O (right 
column), are in good agreement with Ref. [27] . 


4 Charm-quark fragmentation into hadrons in e^p 
collisions 

The charm-hadron-production cross-sections at HERA were measured in a restricted hdu- 
cial phase space. The extraction of the charm-quark fragmentation fractions requires a 
special treatment, as described in detail in App.jB] The approach followed in this analysis 
is similar to the one originally used by the ZEUS collaboration m- 

4.1 Charm-quark fragmentation fractions from measurements 
in DIS 

Charm-quark fragmentation fractions in DIS in e^p collisions are calculated from ZEUS 
and HI measurements given in Tab. [6l 

For the calculation of charm-quark fragmentation fractions a £t procedure is used as 
it is described in Sec. |2J The free parameters in the ht are the charm fragmentation 
fractions and pairs of variables a(c)j|j=i ...3 and for each set of measurement. 

Here, a(c)j is the total charm cross-section in e^p , while Ki is the kinematic factor for 
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Decay 

ZEUS [16] 

ZEUS [IT] 

HI [18] 

CT, nb 

CT, nb 

(T, nb 

D° ^ K+-K- 

7.34 ±0.36lp p± 0.13 


6.53 ± 0.49 

D+ K + TX+-K+ 

2.80 ±0.30t°d|± 0.10 


2.16±0.19l°|® 

D+ K°-k+ 

D+ <f>{K+K-)-K+ 

1 27 + n 1 fi+OH +0.19 
l.Zt m U--LD_o.06-0.15 

25.7 ±4.lt|| ±0.8 

1.67 ±0.41^°;®^ 

D*+ D°{K-tt+)tt+ 

3.14± 0.12to;i5 ±0.06 


2.90 ± 0.2oIo 44 

A.TT+ 


14.9 ±4.9l2.6 ±3.9 


A+ ^ ROp 


14.0 ±5.8t|;|± 3.7 


(combination) 

K-7T+,no 

1.78 ± 0.08lo lo ± 0.03 

14.7±3.8t|2 ±3-9 



Table 6: Measurements of charm-hadron-production cross-sections in DIS in 
e'^p collisions. The numbers are given as in the original publications. The first 
uncertainty is statistical, the second is systematical and the third one corresponds 
to the branching ratio. 

decays from higher states (see App. [B]). The parameter k is hxed to one for the low-p^ 
measurements in Ref. mi since the whole pr kinematic space was covered. The present 
beauty contributions in Ref. [IB] have small impact on the result and are neglected. The 
sum of charm fragmentation fractions Sep dis is constrained to unity. The results of the 
averaging procedure are given in Tab. [T] 



Constrained S 

f(c^D*+) 

0.2372 ±0.0173 

f{c^D+) 

0.2170 ±0.0203 

fic^DO) 

0.6272 ± 0.0287 

f{c ^ Dt) 

0.0945 ± 0.0124 


0.0540 ±0.0195 


1.7 

^dof 

3 

s 

1.0000 ±0.0004 

^u/d 

1.2361 ±0.1331 

pd 

0.6282 ± 0.0440 

Is 

0.2240 ± 0.0320 


Table 7: Average of charm-guark fragmentation fractions in e^p collisions in DIS. 
The guantities S, Ru/d, Py 7s recalculated from the fit results taking into 
account correlation of fit parameters. The value of minimised and the number 
degrees of freedom of the fit Udof are given as well. 


The obtained fragmentation fractions are in agreement with those obtained in the original 
publications [MIITB] . The uncertainties of the obtained results are somewhat larger because 
this analysis, contrary to those studies, relies only on DIS results, whereas the HERA DIS 
papers [IHlilB] used fragmentation fractions into A+ measured at e+e" colliders. 
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4.2 Charm-quark fragmentation fractions from measurements 
in PHP 

Charm-quark fragmentation fractions in PHP in e^p collisions were calculated from mea¬ 
surements of ZEUS collaboration and given in Tab. 13 


Decay 

ZEUS[T4] 

(T, nb 

ZEUS[I5] 
f{c ^ H) 

D+ K-TT+1T+ 

5 Q 7 _|_ n Of?+0.44+0.34 
O.Uf It U-'JD_o.23-0.30 

0.234 ± 0.0061!;“^ 

D° K+TT-,no D*+ 

8 4Q + 0 44+0-47-I-0.20 


D° K+TT-,with D*+ 

2.65±0.08lp;iJ ±0.06 

0.588 ±0.017l!;!“l!'°l! 

D+ 4>(K+K-)it+ 

2.37 ± 0.20 ± 0.2o1:o 45 

0.088 ±0.0061°;°°?!°;°°® 

A+ pK-TT+ 

■3 crq _(_ n f?/^ + 0.54-|-1.15 

O.oy It lJ-OD_o.66_o.70 

f) nyg _(_ n ni o+0.005+0.024 
U.U/y It D.Uid_o. 009-0.014 

D*+ D°{K--k+)-k+ 


It D.UUO_o 004-0.007 

add£)*+ £)0 

1.05 ±0.071“;°! ±0.03 


km£)*+ ^ £,0 

4 97 + 0 14+0-23-1-0.13 
=>= D.i^_o.i8_o.i2 



Table 8: Measurements of charm-hadron-production cross-sections and fragmen¬ 
tation fractions in photoproduction in e^p collisions. The numbers are given as in 
the original publications. The first uncertainty is statistical, the second is system¬ 
atic and the third one corresponds to the branching ratio. 

For the update of the latest ZEUS measurement jlS] to the decay branching ratios from 
Tab.m the measured fragmentation fractions are hrst transformed into total charm-hadron 
cross-sections according to the formulae in App. |B]and only then used in the calculations. 
In this procedure, the kinematic factor for decays from higher states, k, is set to 1, since 



Constrained S 

fic^D*+) 

0.2345 ± 0.0081 

f{c^D+) 

0.2341 ± 0.0093 

f{c^D°) 

0.5991 ±0.0126 

f(c^D+) 

0.0901 ± 0.0062 

/(c —> A^) 

0.0675 ± 0.0106 


5.2 

^dof 

4 

s 

1.0000 ±0.0005 

^u/d 

1.1209 ±0.0545 

jDd 

0.5970 ±0.0181 

'Is 

0.2164 ±0.0162 


Table 9: Average of charm-guark fragmentation fractions in hadrons in e^p 
collisions in photoproduction. The guantities S, Ru/d, Py Is are recalculated 
from the fit results taking into account correlation of fit parameters. The value of 
minimised and the number degrees of freedom of the fit Udoj are given as well. 
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the total phase space is considered from the fragmentation fraction dehnition, and the 
cr(c) value cancels out in the procedure. 

For the calculation of charm-quark fragmentation fractions, a £t procedure is used as 
it is described in Sec. [2J The free parameters in the fit are the charm fragmentation 
fractions and pairs of variables cr(c)j|j=i ,,,2 and Ki\i=i ,,,2 for each set of measurement. Here, 
cr(c)j is the total charm cross-section in e^p, while k* is the kinematic factor for decays 
from higher states (see App. [B]). The sum of charm fragmentation fractions S^p php is 
constrained to unity. The results of the averaging procedure are given in Tab. [9l The 
obtained fragmentation fractions are in agreement with those obtained in the original 
publications [iiliig. 


5 Charm-quark fragmentation into hadrons in pp co¬ 
llisions 

The ALICE experiment measured hducial cross-sections of Df [22] and differential pt 
cross-sections of D^, and D*~^ mesons [2T1I23] at ^/s = 2.76 TeV and ^/s = 7TeV. With 
an integration of the differential cross-sections of D^, and D*~^ from Ref. [21] and 
from Ref. [23], a coherent set of measurements in the kinematic range 2 < pp < 12 GeV, 
IpI < 0.5 has been constructed (see Tab. [10]) for the -y/s = 2.76 TeV and y/s = 7TeV. 
The LHCb experiment provided measurements of charm-hadron cross-sections at y/s = 
7 TeV [IH] and at y/s = 13 TeV [20]. The ATLAS experiment recently measured the 
production cross-sections of D*^, and D'^ mesons at y/s = 7 TeV [23] in the kinematic 
range 3.5 < pp < 20 GeV, |p| < 2.1. 

The measurements together with the correlation matrix for LHCb y/s = 7 TeV are given 
in Tab. flOl 

For the calculation of charm-quark fragmentation fractions a ht procedure is used as it is 
described in Sec. |2J The free parameters in the fit are: the charm fragmentation fractions, 
the fiducial cross-sections for LHCb, ALICE and ATLAS measurements in kinematic 
regions given in Tab. [inland corresponding k parameters. 

The constraint Spp — 1 = 0 is imposed. As the Refs. [2I1E3] do not provide detailed decom¬ 
position of the systematic uncertainties, for every bin all systematic uncertainties were 
conservatively assumed to be fully correlated. For all of these measurements we assume 
the statistical and systematic uncertainties uncorrelated and luminosity uncertainties - 
fully correlated within a set of measurements at a given value of y/s. 

A set of orthogonal fully correlated uncertainties was obtained from the covariance ma¬ 
trix of the y/s = 7 TeV LHCb measurements with an eigenvector decomposition. The 
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LHGb [19] 

LHGb [20] 

ALICE [23] 

ALICE [2T] 

ATLAS [24] 


7 TeV 

13 TeV 

2.76 TeV 

7 TeV 

7 TeV 

Decay 

PTe[0,8] 2/e [2,4.5] 

Pt6 [1,8] 

PTe]2,i2] 

PTe]2,12] 

PT 6 [3.5, 20] 



Corr.(%) 

ye [2,4.5] 

\y\ < 0.5 

\y\ < 0.5 

[r,] <2.1 


(j, /ib 

D° D+ D-+D+ 

CT, fih 

(T, fih 

(T, /ib 

cr, fib 

D° Kit 

1661 ± 129 


2460 ± 3 ± 130 

110 ± 16^34 

231 ± 12t|g 


D+ K-TT+1T+ 

645 ± 74 

76 

1000 ± 3 ± 110 

47± 9ll§ 

81 ± 

328 ± 16 ± 27 

D*+ D°{K-tt+)it+ 

677 ± 83 

77 73 

460 ± 13 ± 100 

59 ± lihll 

104 ± 6tll 

331 ±18 ±28 

Dt <j){K+K-) it+ 

197 ± 31 

55 52 53 

880 ± 5 ± 140 


53± 12l“ [22] 

160 ±31 ±17 

AJ —>■ pK~'K'^ 

233 ± 77 

26 25 25 18 








PT 6 [0,8] 




D° Kit 



3370 ± 4 ± 200 




D+ —>■ K~tt+tt+ 



1290 ±8± 190 





Table 10: Measurements of charm-hadron-production cross-sections in pp colli¬ 
sions. For the measurements from Ref. m the total uncertainty is given. For the 
rest of measurements the first uncertainty is statistical, the second is systematic. 

Pt is gived in GeV. 

obtained uncertainties are later treated in the same way as other correlated sources in 
the combination. The Ref. [20] does not contain the correlation matrix for the mea¬ 
surements at ^/s = 13 TeV, therefore simplihed correlations between measurements were 
calculated as follows. All of the measurements include 3.9% fully correlated uncertainty 
related to luminosity included in the Ref. ESI to the systematic uncertainty. The sys¬ 
tematic uncertainties also include the uncertainties on the branching ratios, which were 
treated correlated with other branching-ratio uncertainties. The remaining systematic un¬ 
certainty were treated as fully uncorrelated for different measurements with the same pt 
cuts and fully correlated for the same measurements with different pt cuts. The statistical 
uncertainties of cT(Zl°)py<8GeV, cr(-D''")pT<8GeV were split in two parts, which correspond 
to pr < 1 GeV and > 1 GeV regions. The later were considered fully correlated to the 
statistical uncertainties of the cr(Zi)°)i<pj,<8GeV and cr(Zi)’^)i<pj,<8GeV measurements. 

For the LHGb and ATLAS measurements of the Df —)■ K^K~'k^ in the limited mass 
windows M[K^K~) the following approach was used. The branching ratios were obtained 
from the integrals over the M{K^K~) line shape that was parametrised as in Ref. [31| with 
the total Df —)■ signal normalised to B{Df —)■ 7r+) = 5.45 ± 0.17% [25] . 

The results are given in Tab. [TJ 

The results of the £t are reported in Tab. [TT] In addition to the values of the frag¬ 
mentation fractions, the fit delivers the inclusive charm-production cross sections in the 
corresponding hducial regions, which have particular interest. Therefore, the values of 
these cross-sections obtained in the global combination with better precision are discussed 
below. 
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Constrained S 

/(c^D*+) 

0.2337 ±0.0175 

f{c^D+) 

0.2274 ±0.0128 

f{c^D°) 

0.6176 ±0.0160 

f(c^D+) 

0.0824 ± 0.0084 

/(c —)■ A^) 

0.0639 ±0.0122 


6.9 

^dof 

7 

s 

1.0000 ±0.0005 

^■u/d 

1.1913 ±0.1012 

r>d 

0.6059 ± 0.0306 

7s 

0.1951 ± 0.0216 


Table 11: Average of charm-quark fragmentation fractions in pp collisions. The 
quantities S, Ru/d, Ry 7* recalculated from the fit results taking into account 
correlation of fit parameters. The value of minimised and the number degrees 
of freedom of the fit Udof are given as well. 

6 Selection of measurements for the extraction of 
fragmentation fractions 


The selection of the measurements for the extraction of fragmentation fractions was done 
according a set of criteria explained below. 

First, the selection is limited to the measurements obtained in the collisions of high energy 
particle beams as it assures an absence of possible matter effects and the charm quark 
production mechanism in these environments is well understood. The measurements of 
charm-hadron production in proton-meson, proton-nucleon and nucleon-nucleon collisi¬ 
ons [381H3] were omitted as those provide results in very specihc production environment 
and energy ranges which cannot be easily compared to the results in other experiments. 

The second criteria of the selection is the precision of the measured quantities: the mea¬ 
surements in collisions with ^/s = 12 — 90 GeV from MARK-II [H], HRS |15] - H9] . 
TPC [SU], TASSO [SI1IS2], JADE [S31IS1], VENUS [SS] and some other experiments have 
very limited precision and are not used for the global combination. 

The third criterion of the selection is the availability of sufficient measurements in the 
given physical environment needed for the extraction procedure. Several results on charm 
production in collisions (e.g. Ref. [SB]) and pp collisions (e.g. Refs. [121I131[S7]) do not 
contain enough simultaneous measurements of hadron production and, therefore, cannot 
be treated independently and/or constrain the results of the combination. 


15 


















7 The global combination 


To check the consistency of the data from different production regimes and also to extract 
the charm-quark fragmentation fractions with high precision, all input measurements 
introduced in the previous sections are used together to produce a global combination. 
As discussed in Sec. 13.11 the A+ measurement by the BABAR experiment [S] is not 
included while obtaining the combined result. The free parameters of the £t are the 
charm-quark fragmentation fractions and pairs of variables cT(c)j|j=i ...5 and Ki|i=i ...5 for 
three DIS and two PHP sets of measurements, a(e~^e~ —)■ c) at y/s = 10.5 GeV, and 
the hducial charm-quark cross-section and k parameters in pp collisions, at y/s = 7 TeV 
and y/s = 13 TeV, corresponding to the phase space of the measurements. The constraint 
on the sum of the cross-sections of the weakly decaying charm states, S, is imposed in 
the combination, i.e. the prediction for the total charm cross-sections in e’''e“ collisions 
is not used, in order to minimise model dependence of the averaging procedure. The 
result of averaging and pp data, with the constraint S' = 1 is presented in 

the middle column of Tab. [T2] and is shown in Fig. [T] The correlations between the 
fitted parameters are given in Tab. [I3l The input data are in very good agreement with 
X^/^doi = 65.6/64. The result of the combination has signihcantly reduced uncertainties 
compared to individual measurements. 



Constrained S 

Constrained S, 
fixed cr(e“^e“ c), 

f{c^D*+) 

0.2429 ± 0.0049 

0.2386 ± 0.0046 


0.2306 ± 0.0315 

0.2250 ± 0.0299 


0.0548 ± 0.0074 

0.0537 ± 0.0072 

fic^D+) 

0.2404 ± 0.0067 

0.2439 ± 0.0067 

fic^DO) 

0.6086 ± 0.0076 

0.6141 ± 0.0073 

nc^D+) 

0.0802 ± 0.0040 

0.0797 ± 0.0040 


0.0623 ± 0.0041 

0.0549 ± 0.0026 


65.6 

87.1 

R'doi 

64 

67 

s 

1.0000 ± 0.0005 

1.0000 ± 0.0004 

^u/d 

1.0971 ± 0.0354 

1.1164 ± 0.0354 

ryd 

0.5578 ± 0.0375 

0.5403 ± 0.0355 

Is 

0.1890 ± 0.0103 

0.1859 ±0.0101 

A 

0.2314 ±0.0347 

0.2316 ±0.0346 


Table 12: Average of charm-quark fragmentation fractions in hadrons. The 
quantities S, Ru/d, Py cmdjs are recalculated from the fit results taking into account 
correlation of fit parameters. The value of minimised and the number degrees 
of freedom of the fit Udof are given as well. 
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Combination 
■ e"^e' Z decays 
A e±p PHP 


• e^e' B-factories 

♦ e±p DIS 
T pp LHC 



Figure 1: The values of charm-quark fragmentation fractions, /(c —)■ H), in 
different experiments with the S constraint. The global combination with the S 
constraint is shown with the shaded band. Averages of included data in different 
production regimes are shown with various full symbols. 

As an alternative, the combination is also performed using both the constraint on S 
as well as theoretical predictions of charm production in e’''e“ collisions and Z decays, 
i.e. a{e~^e~ —)■ c) at ^/s = 10.5 GeV and This approach also allows to include the 
precise BABAR measurement of A+ production [S] using the calculation as described in 
App-O which signihcantly affects the averaged value of /(A+). The result of the averaging 
procedure with this approach is given in the right column of Tab. [T2] for completeness. 
The result is more model dependent than the default combination, but has a higher 
precision. At the same time, the result visibly differs from the result of the default 
procedure. This may partially be traced to the value Sz = 0.9292 ± 0.0261 for the 
accurate LEP measurements, which differs markedly from 1 (see Sec. 13.21) . This difference 
is also reflected in the larger y^/udof value compared to the default global combination. 
The difference in the f{c —?■ A+) precision is to a large extent due to inclusion of the 
precise BABAR data [S]. 
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D*+ 

to 

o 

Dt+ 

D+ 

D° 

D+ 

At 

D*+ 

1.00 

-0.02 

-0.02 

-0.08 

0.19 

-0.07 

-0.12 

D*o 

-0.02 

1.00 

0.02 

-0.07 

0.07 

0.01 

-0.01 

Dt+ 

-0.02 

0.02 

1.00 

-0.05 

-0.07 

0.23 

-0.01 

D+ 

-0.08 

-0.07 

-0.05 

1.00 

-0.66 

-0.19 

-0.19 

D° 

0.19 

0.07 

-0.07 

-0.66 

1.00 

-0.32 

-0.41 

D+ 

-0.07 

0.01 

0.23 

-0.19 

-0.32 

1.00 

-0.07 

A+ 

-0.12 

-0.01 

-0.01 

-0.19 

-0.41 

-0.07 

1.00 


Table 13: Correlation of charm-quark fragmentation fractions from the fit with 
constrained S. 



A 

O 


Combination 
e*e' Z decays 
e±p PHP 
pp CDF 


• e*e' B-factories 

♦ e±p DIS 

T pp LHC 

□ Pv average by A. David 



Figure 2: The values of Ru/d, Py Is in different experiments with the S 
constraint. The global combination with the S constraint is shown with the shaded 
band. Combinations of included data in different production regimes are shown with 
various full symbols. Data that were not included in the combination are 

shown with open symbols. Note, that the latter are quoted from the original papers, 
i.e. without correction to the up-to-date branching ratios and with no branching 
ratio uncertainty, if not given in the source. 
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The extracted Ru Id, Py and 7 <j factors are provided in Tab. [12] and shown in Fig. [2] The 
combined data are also compared to recent measurements [53,[58] that were not included 
in the combination. In particular, Ru/d = 1-097 ± 0.035 is in fair agreement with the 
isospin invariance hypothesis Ru/d = 1 within 2.7 standard deviations. The values of 
the a{pp —)■ c) cross-sections, obtained in the global fit (see Tab. [Tl|) are consistent with 
those obtained in the original analysis, but have significantly reduced uncertainties. The 


Vs, 

TeV 

PT 

range, 

GeV 

y or T} 

range 

Fit result 

cr(pp c), 

fib 

Original 

cr{pp c), 

fib 

7 

[0,8] 

P 6 [2,4.5] 

2689 ± 203 

2838 ± 268 [l9] 

13 

[1,8] 

P 6 [2,4.5] 

4174 ± 339 

4300 ± 356 [20] 

13 

[0,8] 

P 6 [2,4.5] 

5269 ± 293 

5880 ± 482 [20] 

2.76 

[2,12] 

1p1 <0.5 

229 ± 67 


7 

[2,12] 

\y\ < 0.5 

434 ± 84 


7 

[3.5,20] 

Id <2.1 

1400 ± 141 



Table 14: The values of the inclusive fiducial charm quark production cross- 
section, a{pp c), from the original publications and obtained in the global fit. The 
statistical, systematic and fragmentation uncertainties of the values from Refs. fIR 
were added in quadrature. 


consistent treatment of the LHCb and ALICE measurements in the combination procedure 
allows unbiased calculation of the ratio of the inclusive fiducial charm-quark production 
cross-sections: 


-R 7 / 2.76 


a{pp —)■ c)7TeV,2<pT<12GeV,|i/|<0.5 
a{pp —>• c)2.76TeV,2<PT<12GeV,|y|<0.5 


1.89 ± 0.66 


and 

a{pp c)7TeV,pT<8GeV,2<P<4.5 

The Rt 12.116 is compatible with the predictions in Ref. [59]. The R 13/7 value is visibly 
higher than the theoretical prediction Ri 3 / 7 (theory) = 1.39 ^q;29 


8 Excited states 

In addition to the average fragmentation fractions for the ground, L = 0, states, some 
fragmentation fractions for the excited, L = 1 charm hadrons are calculated. 

The measurements used for the calculations are shown in Tab. HSJ The unpublished 
measurement of f{c —)■ Hj"!) from Ref. [M] was not used. The fragmentation fractions 
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were not updated to the most recent branching ratios, as the difference between the used 
branching ratios and the newest is negligible in comparison to statistical and systematical 
uncertainties of the measurements, and is well below the given numerical precision of the 
individual measurements in Tab. m 



ZEUS 

OPAL 

ALEPH 

Particle 

/(c ^ H) 

f{c ^ H) 

f{c ^ H) 


(10-2) 

(10-2) 

(10-2) 


2.9±0.5tp 5 E2] 

3.5 ± 0.4to 6 [33] 

2.1 ±0.7 ±0.3 1631 


ry*Q 

^2 

3.9±0.9lp g [62] 
3.8±0.7t°;g [33] 

5.2 ± 2.2 ± 1.3 l63l 


Df 

4.6± l.Slg 3 [62] 



p* + 

-^2 

3.2±0.8t°;2 [62] 



Dti 

l.ll±0.16l° °® [33] 

1.6 ± 0.4 ± 0.3 1631 

0.94 ± 0.22 ± 0.07 [64] 


Table 15: Comparison of fragmentation-fraction-results of measurements of ex¬ 
cited charm mesons. 

The averages are calculated with an assumption of fully uncorrelated statistical and sys¬ 
tematical uncertainties. The results of the averaging procedure are given in Tab. [161 The 



Average (10 2) 

f(.c^D+) 

4.60t?;S® 

f{c ^ D*+) 

o 9|-»+0.94 

g2 

f{c ^ DO) 

2.97 ±0.38 


3.94 ±0.68 


1.09 ±0.14 

7s 1 

Qo -^+8.0 


Table 16: Average of charm-guark fragmentation fractions in excited charm 
mesons. The guantity is calculated from the averaging results without taking 
into account correlations. 


strangeness-suppression factor for L = 1, J = charm mesons is calculated neglecting 
11(2430)*^ contribution and assuming Df is state: 


7.1 


2/(c ^ Dtd 

f(c^D<l) + !(c^Dty 


9 Summary 

A summary of measurements of the fragmentation of charm quarks into a specific charm 
hadron is given. The analysis includes data collected in photoproduction and deep inelas¬ 
tic scattering in collisions and well as e’''e“ and pp data. Measurements in different 
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production regimes agree within uncertainties, supporting the hypothesis that fragmen¬ 
tation proceeds independent of the specific production process. Averages of the fragmen¬ 
tation fractions are presented. The global average has signihcantly reduced uncertainties 
compared to individual measurements. In addition, the hypothesis that the sum of frag¬ 
mentation fractions of all known weakly decaying charm hadrons is equal to unity is 
checked to hold within 3 standard deviations using the e'''e“ data. 
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Appendix A Predictions for charm production at B- 

factories 

The total cross-section of quark q production in e’''e“ collisions at energies 2M{q) <C 
^/s M{Z) can be given as 

a{e^e~ q) = 2a{e'^e~ —)■ ^ (V 

colours 

where Vg is the vector electromagnetic coupling of the quark q (i.e. charge), rq(s) are the 
correction coefficients with higher order QCD corrections and 

a(e~^e~ — )■ = 4a^(s)7r/3s ( 8 ) 

is the total cross-section of massless charged lepton pair production. In this work, the 
calculations of the rg(s) were done according to Ref. [31] at the reference energy of 
= 10.5 GeV and assuming the c quark is the heavy one. The constants used for the cal¬ 
culations in Eq. ([7]) and Eq. (jH]) are the strong coupling as(\/s = 10.5 GeV) = 0.172 [21], 
the MS charm-quark mass mc{\fs = 10.5 GeV) = 0.74 GeV [3T] and the electromagnetic 
coupling a{^/s = 10.5 GeV) = 1/132.0 (calculated according to [251ISS] as implemented 
in IS71). The uncertainties on the given values are negligible. The result of the calcula¬ 
tions is a{e^e~ c^\fs = 10.5 GeV) = 2399.23 nb. To verify the calculations, Eq. ([7]) 
can be rewritten as 

a{e^e~ c) = 2a{e^e~ —)■ /^/“)Rc-Rhad, 
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where the quantities 


^ cr(e’''e —)■ hadrons) 
^ a(e+e- ^ /+/") 


njr, = 3.5239 

u,d,s^c colours 


and 


Rc 


a{e'^e cc) 
(T(e+e“ —hadrons) 


E '^c^c 

. = 0.4012 

E E 

u^d^s,c colours 


can be compared with the existing measurements and predictions. It was found that 
i?had is in agreement with the direct measurement from CLEO below y/s = 10.56 GeV 
.Rhad.cLEO = 3.591 ± 0.003 ± 0.067 ± 0.049 [68] and Rc is in agreement with the CLEO 
Monte-Carlo based estimation i?c,CLEO = 0.37 ± 0.05 [T]. For all the theoretically calcu¬ 
lated values, the uncertainties of calculations are negligible. 


Appendix B Extraction of the charm-quark fragmen¬ 
tation fractions from the measurements 
in the restricted phase space 


Often the production cross-sections of charm hadrons are measured in some restricted 
(e.g. in transverse momentum and pseudorapidity) kinematical region v and cannot be 
used directly (i.e. without extrapolation to the full kinematical space) in the Eq. (|2|). To 
avoid the extrapolation and obtain unbiased charm-quark fragmentation fractions, the 
following approach is used. 


The full cross-section of a charm hadron, a{H)^^ can be split into cross-section of direct 
production cT(i7)dir,eij the contribution from the decays of heavier charm states H*, 

^ (77) decays,Ef 


— 0'{H)dir,ev + 0'(77)decays,ei; — 0'{H)dir,ev + ^ ^ 

allH* 

where kH*^H < 1 is a fraction of H* —)■ H decays with H in v. The lack of experimental 
data in this analysis allows to consider only the heavier charm states that are giving the 
largest contribution to the total cross-sections. For this reason all and produced not 
in D*^ and D*^ decays are considered as produced directly. Because of similar kinematics 
oi D* ^ D decays it is also assumed that k£,*^£) = k. The measurements on heavier 
charm baryons and charm strange mesons are absent, but all of those decay dominantly 
to A+ and 77+ so we treat all A+ and 77+ as produced directly. Within a common phase 
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space, the impact of the different masses of the hadrons on the fragmentation process is 
neglected. 

With the assumptions above, we have the equations: 


^ a(D+)6, 

= a(D+)dir,e. + A:a(D*+)dirS(D*+ ^ D+) 

ap°)e. 

= + fca(D*0)^i^D(D*o ^ DO) 

+A:a(D*+)diri3(D*+ ^ DO) 


C^i^s )dir,ew 

c^(A+)e. 

‘^(Ac )dir,ew 

a(D*+)6, 

= a(D*+)dir,e,; 

^ a(D*0),, 

= a(D*0)di^g^. 


Assuming cr(if)dir,eij = c’'(c)et;/(c —t hf)dir and introducing k = ka{c) ja{c)^.u we have: 


= 


In the full kinematical space: 

' f{c^D+) 

f{c ^ D^) 


fic ^ Dt) 
fic ^ A+) 
fic ^ D*+) 
. fic ^ D*^) 


cr(c)6^(/(c /^+)dir 

+Kfic ^ D*+U,BiD*+ ^ D+)) 

cr(c)6„(/(c -)■ £>°)dir 

+Kfic D*^)dnBiD*^ D^) 

+Kfic ^ D*+U,BiD*+ ^ DO)) 

cr(c)6^/(c D + )dir 

cr(c)6^/(c A+)dir 

Cr(c)6^/(C D*+)dir 
a(c)6^/(c D*0)dir. 


— /(c —!■ D'^)dir 

+/(c ^ D*+U,BiD*+ ^ D+) 

= /(c —)■ D0)dir 

+/(c ^ D*0)diri3(D*0 ^ DO) 
+/(c ^ D*+)dirS(D*+ ^ DO) 
= /(c—tD + )dir 
= /(C—)-A+)dir 

= /(c^D*+)dir 

= fic^D*^U,. 


In general, to solve the system the measurements of D*o production are needed. However, 
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these can be avoided with an assnmption of isospin invariance: 


/(c —)■ -D''')dir _ /(c —)■ -D*’'')dir 
/(c £>°)dir /(c £>*°)dir ’ 

The last two systems are the working equations for the calculation of the charm fragmen¬ 
tation fractions from the cross-section measurements in the restricted phase space. 
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